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Motivated by recent claims of lines in the Fermi gamma-ray spectrum, we critically exam- 
ine means of enhancing neutralino annihilation into neutral gauge bosons. The signal can 
be boosted while remaining consistent with continuum photon constraints if a new singlet- 
like pseudoscalar is present. We consider singlet extensions of the MSSM, focusing on the 
NMSSM, where a 'well-tempered' neutralino can explain the lines while remaining consis- 
tent with current constraints. We adopt a complementary numerical and analytic approach 
throughout in order to gain intuition for the underlying physics. The scenario requires a 
rich spectrum of light neutralinos and charginos leading to characteristic phenomenological 
t-H signatures at the LHC whose properties we explore. Future direct detection prospects are 

excellent, with sizeable spin-dependent and spin-independent cross-sections. 
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INTRODUCTION 



There is robust evidence on astrophysical and cosmological scales for the presence of particle 
dark matter (DM) in our universe PQ ]2j. A particularly well studied candidate is the neutralino 
and determining the properties of this particle is the subject of an intense experimental effort 
at colliders, direct detection and indirect detection experiments. Indirect detection experiments 
search for particles created from the decay or annihilation of DM particles from regions where the 
density of DM is expected to be high, such as the centre of our galaxy. A particularly striking 
signal is a monochromatic 7-ray line, arising when the DM annihilates into a two-body final state 
containing a photon [3]. 

Recently, numerous studies of the publicly available Fermi-LAT [4] data found a sharp fea- 
ture in the 7-ray spectrum at ~ 130 GeV coming from the vicinity of the galactic centre [5HS]. 
Interpreting the feature as a monochromatic line arising from DM annihilation into two pho- 
tons, Weniger [6] found that DM with mass 129.8 ± 2.4±i3 GeV and annihilation cross-section 
(av) rr = (1.27 ±0.321^1) x 1(T 27 cmV 1 fits the signal well (see also [7-14]). Intriguingly, 
Refs. [5J [T3] note that two lines, one at ~ 130 GeV and a weaker one at ~ 111 GeV, provide 
a slightly better fit to the Fermi data. Such a pair of lines can be naturally explained by a DM 
particle of mass ~ 130 GeV annihilating into 77 and with a relative annihilation cross-section 
(crv)- f z / (w)^ = 0.66iQ4g [16] . In this paper, we address under what conditions the neutralino 
can fit these observations and the implications for the neutralino sector. 

Doubts about a DM origin of these features have been raised in |17H19j . although as yet there 
is no compelling astrophysical process that can explain the Fermi features. Searches for systematic 
effects associated with the Fermi-LAT instrument have been performed in [20-22J and a small 
excess of photons from the Earth's limb at ~ 130 GeV for photons within a limited range of 
detector incidence angles was found. However, such an effect cannot account for all of the signal 
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from the galactic centre [2H 122] . Finally, the annihilation cross-section required to explain the 
feature is in mild tension with the upper limit (<tw) 77 < 1.0 x 10~ 27 cm 3 s _1 set by the Fermi-LAT 
Collaboration from the 2 year dataset [23] . Furthermore, preliminary analysis of the 4 year dataset 
finds a feature at ~ 135 GeV but with a lower significance compared with the analyses quoted 
above |24j . The fate of the feature will ultimately be resolved as more data is collected. For the 
purposes of this paper, we assume that the features observed by Fermi are due to photons arising 
from neutralino annihilation and study the associated phenomenological implications. 1 

Interpreting the feature as a DM signature has its own challenges: Ref. [8] find that the feature 
is offset from the galactic centre by 1.5°, although this may arise due to an interplay between 
the galactic baryons and DM [26]. A more serious issue is explaining the relatively large size of 
the annihilation cross-section into photons that is required while remaining consistent with the 
continuum flux of photons arising from annihilation of the DM into W and Z bosons and Standard 
Model fermions |25t [2"7rl30| . As the DM is a neutral particle, its coupling to photons generally 
arises at loop-level and is suppressed relative to the coupling to W and Z bosons and Standard 
Model fermions, which can occur at tree-level. Ref. |29| quantified this constraint through the 
quantity 

R th = - } av) ^ - , (1) 
2{av) 11 + {(tv)jz 

where {crv) arm is the total annihilation cross-section and (av)~ i ~ it ~ i z are the annihilation 
cross-sections into 77 and 7Z respectively. When the total annihilation cross-section is 
(<7t>)ann ~ 3 x 10 -26 cm 3 s _1 (required to achieve the correct thermal relic density) and when 
(<Tf) 77 ~ 1.2 x 10~ 27 cm 3 s _1 and {crv) 1 z / ' (0^)77 = 0.7 (to explain the lines in the Fermi-LAT 
data), then i? th ~ 9, compatible with all continuum constraints in |29j . Thus, if the thermal 
relic density and line-strength for a given DM candidate are correct, it should be compatible with 
continuum photon constraints. 2 

Many papers have been written with methods and models that explain the monochromatic fea- 
tures while remaining consistent with known constraints, see e.g. |27U33fJ53] . Within the framework 
of the simplest supersymmetric extension of the Standard Model, the MSSM, the large annihila- 
tion cross-section into photons (crv) 77 can be achieved |54J. However, this requires a wino- or 
higgsino-like neutralino and such scenarios yield a negligible thermal relic density and also give 
rise to a large continuum flux [55], which rule them out [23 [281 Said another way, if the 
MSSM neutralino achieves the correct thermal relic density then (cn>) 77 is orders of magnitude too 
low to explain the Fermi feature. 3 Therefore, we see what is required: we need to keep the total 
annihilation cross-section at the level required to obtain the correct relic density while boosting 
the cross-section into photons. This can be achieved by exploiting a resonance with a pseudoscalar 
that couples primarily to photons, as suggested in (27J E7J [58] • Although the MSSM does have 
a pseudoscalar Higgs A, this approach fails because A does not couple primarily to photons; it 
couples to charged states // at tree-level, while the coupling to photons is generated at loop level. 
However, as discussed in |57| |59"1 160j . within extensions of the MSSM that include extra singlet 
states, the tree-level coupling to // can be reduced while maintaining the coupling to photons. 



1 Decaying DM is disfavoured as it requires an enhancement of the DM density near the galactic centre |25| . 

2 Similar conclusions can also be drawn from the constraints from pp and synchrotron radiation. Again, if the 
thermal relic density is correct, the DM candidate is compatible with these constraints, see e.g. |25l 131] and [32] . 

3 In fact, even in the MSSM the 130 GeV line can be explained by internal bremsstrahlung from a bino-like neutralino 
annihilating to light leptons .16, 56^ . In order to enhance the number of photons produced, an accidental degeneracy 
between the mass of the exchanged t-channel slepton and the neutralino is required. We will not consider internal 
bremsstrahlung further in this paper. 
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Singlet extensions of the MSSM typically extend the Higgs and neutralino sectors of the MSSM 
and have received much attention as they can help to reduce the electroweak fine-tuning (see 
e.g. [61H66] ) . Indeed, Refs. |59| [60] have shown that the correct annihilation cross-section into 
photons can be achieved in the (Z3-symmetric) NMSSM while achieving the correct thermal relic 
abundance and having a continuum photon flux consistent with observation. In those papers the 
neutralino is mostly bino, with a sub-dominant higgsino component. However, the benchmark point 
given in [!59] predicts a neutralino-nucleon cross-section that is now excluded by the XENON100 
Collaboration's search with 225 days of data [67]. 

In this paper, we reexamine the (Z3-symmetric) NMSSM case in detail and again find that the 
neutralino is dominantly bino- higgsino- like. As well as a numerical implementation, throughout, 
we strive for an analytic understanding of the phenomenology. In Section [TT] we discuss the general 
features that must be present for the neutralino to explain the lines and show that these are present 
within singlet extensions of the MSSM. We then consider the constraints from the requirement of 
obtaining the correct thermal relic density and from direct detection experiments. In particular, 
we describe how this scenario can be made consistent with the XENON100 bound. In Section ITTT1 
we present three benchmark points with distinctive phenomenological signals and consider future 
tests of the model, including signals at the LHC. In Section IV we briefly discuss the neutralino 
phenomenology in a more general singlet extension of the MSSM that does not have a discrete Z3 
symmetry before concluding in Section |V} Finally, details of the numerical procedure we follow to 
calculate (cu) 77 and (av) 1 z as well as various analytic results are presented in appendices. 



II. NEUTRALINO DARK MATTER IN SINGLET EXTENSIONS OF THE MSSM 

The superpotential of the NMSSM, an extension of the MSSM by a gauge singlet chiral superfield 
S and a discrete Z3 symmetry is 4 

VVnMSSM = VlVukawa + \SH U ■ Hd + -S 3 . (2) 

Here, WVukawa is the usual MSSM superpotential containing the Yukawa couplings. The additional 
soft terms associated with the gauge singlet scalar are 

AV soit =m 2 s \S\ 2 + (\A x SH u -H d + ^S 3 + h.c)j . (3) 

In extensions of the MSSM by one singlet superfield, the neutralino mass matrix is diagonalised 
by an orthogonal real matrix N in a similar fashion to the MSSM neutralino mass matrix. The 
resulting eigenvalues are then real, but not necessarily positive [68]. The composition is 

X? = N n B + N 12 W + N 13 H d + N lA H u + N 15 S . (4) 

Following standard convention, we define the bino- fraction as N^, the wino- fraction as N 2 2 , the 
higgsino- fraction as N 2 3 + iV^ 4 and the singlino- fraction as N^ 5 . 

Boosting the 77 rate can now be easily achieved: A pure-singlet pseudoscalar Higgs As has no 
tree-level couplings to Standard Model fermions or gauge bosons. The dominant interaction is with 
neutralinos and charginos through the additional XSH U ■ H d and k5 3 /3 superpotential terms. The 
charginos are necessarily heavier than the lightest neutralino so that, when 2m^o ~ ttia, the (non- 
relativistic) annihilation of neutralinos into photons is resonantly enhanced, through the diagrams 
shown in Fig. [TJ 



4 Our notation is such that H u ■ Hd = H« H d — H^H^- 
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FIG. 1: The left and right panels show the dominant diagrams contributing to (av) 11 and (av) 1 z when 
2m^o ~ rriA s . The s-channel pseudoscalar Higgs As is singlet-like in order to suppress the contribution 
to the continuum flux of photons. Both charginos contribute to the process and in the right panel, the Z 
coupling can mix contributions from different charginos \f an d xf ■ 



In the NMSSM, requiring a pseudoscalar with ra^ s ~ 2m^o means that the singlet component 
of a neutralino with m^o ~ 130 GeV is small: A singlet-like pseudoscalar Higgs As is obtained 
in the NMSSM by requiring that the (1,2) element of the CP-odd mixing matrix vanishes. This 
occurs if A\ ~ 2ks. The resulting mass of the singlet-like pseudoscalar is 

m\ s ~ -3kA k s = -3 - A K fi cS , (5) 

where we have used // e ff = As. If A K ~ fi c Q ~ 0(100 GeV), we require k/X ~ 1 in order that 
rriAa ~ 2m^o. The singlino mass is 

mg = Iks = 2 (^) ji eS > fi cS (6) 

and since this is larger than /x e gf , a neutralino with mass ~ 130 GeV will have a small singlino 
component (in the NMSSM examples we present below it is less than 1%). 



A. The monochromatic line from 77 



The annihilation cross-sections for XX ~ > 77 an d XX ~ > 1% were computed in the MSSM in \<o9\ 
0] and adapted to the NMSSM in |57| . We have implemented a complete numerical calculation 
of these cross-sections with the SloopS code [Ml fTTH73"] , which we describe in Appendix [Aj In this 
section we discuss the process XX ~^ 77 > leaving our discussion of XX ~^ 1% to Section III A In 
order to gain intuition for these processes we also take an analytic approach. We have calculated the 
contribution to the cross-sections from the diagrams shown in Figure [TJ which dominate when As 
is mostly singlet and near the resonance txia ~ ^wi^o. The full expressions are given in Appendix |B| 
We compare the numerical and analytic results in Fig. [2] and find good agreement (this is for the 
'well-tempered' benchmark point, discussed in detail in Section III). To simplify the full expression 
and highlight the parameter dependence of (a"u) 77 , we consider the result assuming that only the 
lightest chargino contributes 
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FIG. 2: The annihilation cross-section (<ru) 77 as a function of the pseudoscalar Higgs mass vtla in the vicinity 
of our 'well-tempered' benchmark point (discussed in more detail later). The solid red line is the result of 
our analytic formula (Eq. (Bl|) including only the diagrams with the pseudoscalar As- The dotted blue 
line shows the result of the full numerical calculation, detailed in Appendix |Aj the agreement is good. The 
width of the resonance at (av)^ ~ 10~ cm 3 s _1 is ~ 1 GeV. 



Here, a is the EM fine structure constant, U and V are the usual matrices that diagonalise the 
chargino mass matrix and m~+ is the chargino mass. In the second line, to further simplify the 
expression, we have assumed that the lightest chargino is degenerate in mass with the neutralino, 
U\2 = V\2 = 1, T^ = and we have denned 5 X = rriA — 2m^o. The numerical prefactor is the 
line-strength found in [6]. 

In order to boost the signal without tuning the neutralino and pseudoscalar masses too much, 
we desire at least some of the following to be true: 

• A large value of A. 

• A large value of n. 

• A neutralino with a large higgsino and/or singlino- fraction. 

• A light higgsino-like chargino with mass m-± ~ m^o. 

In the NMSSM it is usual to keep A < 0.7 in order that it remains perturbative up to the GUT 
scale. However, the singlet contribution to the lightest Higgs mass is proportional to A 2 sin 2 2/3 so a 
larger value of A may reduce the electroweak fine tuning. This scenario in which part of the theory 
becomes non-perturbative has been dubbed A-SUSY [65^171] and counter to naive expectation, it 
has been shown that gauge coupling unification may be improved [75]. We consider points with 
A ~ 0.7 and above. 

The pseudoscalar As couples to the charginos through the \AsH+ HJ term so the signal is 
enhanced for a light higgsino-like chargino. Having a large higgsino-fraction for the neutralino 
helps in this respect as it implies that there should also be a higgsino-like chargino reasonably 
close in mass (in the limit of a pure higgsino-like neutralino, the neutralino and chargino are mass 
degenerate) . Finally, we argued above that the singlino component of the neutralino is small in the 
NMSSM. Therefore, in this case, the higgsino component should be maximised while remaining 
consistent with constraints from the relic density and direct detection experiments. We address 
these points in turn. 
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FIG. 3: The results of a scan showing the values of \Nis\ and | A^i4 1 that give a thermal relic density in the 
range 0.09 < f^DM^ 2 < 0.13 for a bino-higgsino like neutralino with 125 < m^o < 135 GeV. The parameters 
that were scanned over are given in the text. The higgsino-fraction can be large, ranging from ~ 4 — 25%. 



B. Thermal relic density 

In order for the lightest neutralino to be the DM, it should have the correct relic density. We 
assume that the neutralino makes up all of the DM and that it achieves it abundance through the 
usual mechanism of thermal freeze-out in the early universe |76H79| . We have argued for a large 
higgsino-fraction so we first consider the thermal relic density of pure higgsino neutralino DM. To 
a good approximation this is given by [$QJ 



nDMk2 ~ °- 1 (rTev) ' (9) 

where m^o = /i e ff in this case. For a 130 GeV neutralino, the thermal relic density is too low for 
a pure higgsino to saturate the observed value. This is because higgsino neutralinos couple too 
efficiently to gauge bosons so in order to boost the relic density, we need to reduce the efficiency of 
this interaction. This can be achieved most easily by introducing a bino or singlino component. The 
wino component also couples efficiently to gauge bosons so this should be small in order to maximise 
the higgsino part. Therefore, to maximally enhance the (<rw) 77 in the NMSSM, where the singlino 
component is small, the neutralino should be bino-higgsino-like. To gain intuition for how large the 
higgsino-fraction can be, in Fig.[3]we show the values of | iVi3 1 and |iVi4| for which the observed relic 
density satisfies 0.09 < n DM h 2 < 0.13 and 125 < m^o < 135 GeV in the NMSSM. We have scanned 
over the parameters -0.7 < A < 0.5, -1.1 < k < -0.3, 1.8 < tan/3 < 5, -170 < fi c s < -140 GeV, 
-100 < A K < 200 GeV and 125 < M\ < 145 GeV using NMSSMtools 3.2.0 [HTl |82] , which uses 
MicrOMEGAs |83H86j to calculate the relic abundance. 5 We set A\ = 2ns so that the pseudoscalar 
As is singlet-like. This range captures points in the vicinity of our 'well-tempered' benchmark 
point (see Table [I]). The higgsino-fraction can be large, ranging from ~ 4 — 25%. 



We have not included the \X As — > 77 process in our relic density calculation but it yields a subdominant 
contribution to the total annihilation cross-section that sets the relic abundance: (av)-,-, ~ 0.03(<ri;}ann- 
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C. Direct detection 



DM direct detection experiments search for low energy collisions between DM and a target 
nucleus. The low energy effective Lagrangian describing the neutralino-quark interaction is 

£ c ff = aqXiXiM + dqXil^lbXiqi^im , (10) 

where we only consider contributions that are not velocity or momentum-transfer suppressed. 
The first-term gives rises to a spin-independent interaction while the second gives rise to a spin- 
dependent interaction. The spin-independent interaction is mediated by CP-even Higgs exchange. 6 
In singlet extensions of the MSSM, the extra CP-even Higgs boson(s) and the XSH U ■ and 
kS 3 /3 superpotential terms lead to differences from the usual MSSM result for spin-independent 
scattering, as has been discussed in [SBl - fSU] . The extra CP-odd Higgs boson induces a new spin- 
dependent interaction, however, it is momentum-transfer suppressed [90| and negligible compared 
to the dominant interaction mediated by the Z. 

The Lagrangian term a q , responsible for the spin-independent interaction, is proportional to 
the neutralino-neutralino-Higgs coupling gtnxx> which in the NMSSM is given by 

a q oc g hiXX = g(N 12 - tan8 w N^iSuN^ - S i2 N u ) (11) 
- y/2\{SiiN u N 15 + S i2 N 13 N 15 + S i3 N 13 N 14 ) (12) 
+ V2 K S i3 N* 5 . (13) 

The first line reduces to the MSSM result while the second and third lines come from the XSH U ■ Hd 
and kS 3 /3 superpotential terms respectively. The index i runs from one to three over the CP-even 
Higgs bosons, and S is the matrix that diagonalises the Higgs mass-squared matrix. In Appendix [C] 
we provide full details of how the scattering cross-section <7si is calculated from ct q . 

For the case for a mostly bino-higgsino-like neutralino the MSSM- like contribution and \SH U -Hd 
contribution are dominant. We find that there is the possibility of a cancellation in the MSSM- 
like terms between the MSSM-like Higgs contributions when sgn(7Vi3)=sgn(iVi4) |91H93j . Taking 
an explicit example, the 'well-tempered' benchmark point in Table [I] has sgn(iVi3)=sgn(./Vi4) and 
we find Ogj = 1.4 x 10~ 9 pb. Repeating the calculation with the same parameters but with 
iVi4 — > — iVi4, we find Ugj = 13.4 x 10 -9 pb, an order of magnitude larger. The current XENON100 
limit is Ogj = 3.1 x 10~ 9 pb for m^o = 130 GeV [67] 7 so we find that we require sgn(iVi3)=sgn(AT 1 4) 
to be phenomenologically viable. The benchmark point considered in [59J had sgn(Ni3) ^ sgn(Ni4) 
and is now excluded by the 225 live day XENON100 result. 

To understand when sgn(iVi3)=sgn(iVi4), we consider the MSSM tree-level result 

N\ 3 —^es{Mi — m^o) + m| sin /3 cos j3 sin 2 9\y 



Nu mjo (Mi — m^o ) + m| cos 2 f3 sin 2 



in- 



valid in the limit that M 2 is large. This relation is also approximately true in the NMSSM for the 
bino-higgsino case (and when |A|, |k| < 1). In order that sgn(Ni3) = sgn(Ni4), we find 

m 2 7 sin B cos 8 sin 2 9w nn _ , I 5 GeV \ /sin/3cos/3\ 

Aicff < 7T7 — — = 88 GeV — . (15) 

P (M 1 -m x o) [Mi-m x oJ\ 0.23 J v ; 



6 Throughout our discussion we work in the heavy squark limit and ignore the contribution from squark exchange. 

7 This limit is fairly robust against astrophysical uncertainties for DM in this mass range. See e.g [941 195] . 
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For a mixed bino-higgsino-like neutralino, M\ ~ |/x e ff| and M\ > m^o. Thus, we must have fi c s < 0, 
since |^ e ff| < 100 GeV is ruled out by searches from LEP [96J. 

We next consider the spin-dependent interaction, which in the MSSM and its singlet extensions 
is mediated by the Z boson. The Lagrangian term d q oc | iVi 3 1 2 — | A r i4 1 2 , so a spin-dependent 
detection directly probes the higgsino component. In Appendix [Cj we show that the cross-section 
to scatter off a proton is 

cf .4.oxio-% b (H^_!^i!!) 2 . (16) 

The cross-section to scatter off a neutron has the same scaling but the pre-factor is 3.1 x 10~ 4 pb. 

Currently, the best limits on <7p D are from Super-Kamiokande [U7] and the combined IceCube 
and AMANDA-II limits [98J, which have a similar sensitivity at 130 GeV, o| D < 4 x 10~ 4 pb. 
This limit assumes that the capture rate has reached equilibrium and the neutralino annihilates 
solely into W + W~ in the Sun. This is a reasonably good approximation when the higgsino com- 
ponent is large. The expected sensitivity at ~ 130 GeV of the completed IceCube detector 
is o| D « 3 x 10" 5 pb [98] ■ Direct limits can be placed from terrestrial direct detection exper- 
iments but their sensitivity is lower: for instance, COUPP set a limit of <7p D = 6 x 10~ 3 pb at 
130 GeV [99J. The strongest published limit on crjp is from ZEPLIN-III, who find fff <lx 10~ 2 pb 
at 130 GeV [lOOj. XENON100 should improve considerably on this with their 225 days dataset [67]: 
A recent independent analysis |101j obtains a limit of cr,^ D ~ 9 x 10 -4 pb for the same mass. 



III. BENCHMARK POINTS AND FUTURE SIGNATURES 

We now present details of three benchmark points that satisfy all constraints. All points have a 
large enough value of {av ) 77 to explain the Fermi line at ~ 130 GeV, through a resonance with the 
pseudoscalar Higgs As- The important parameters for achieving the correct 7-ray line position and 
strength, as well as the those required to have a singlet-like pseudoscalar Higgs are listed in Table [IJ 
with other more incidental ones detailed in the caption. All points have small values of tan (3 to 
enhance the A 2 tree-level contribution to the Higgs mass; for all points 125 GeV < < 126 GeV. 
We have also chosen all the points to have large values of |A| > 0.7 and based on the results of [75] . 
are unconcerned about A or k becoming non-perturbative between tens of TeV and the GUT scale. 

The first point is a 'well-tempered' point and has the largest higgsino-fraction (19%) of any of 
the points we consider. The large higgsino-fraction results in a large spin-dependent cross-section, 
sitting on the limit of o"p D from IceCube and just below the limit of cr^ D from XENON100. Thus, 
this point should easily be detectable as more data is collected by those experiments. The mass 
difference between the LSP X\ an d NLSP xf ^ s smai b a t ~ 15 GeV. With no intermediate mass 
sleptons, this point is also characterised by three-body decays of X2 anc ^ Xi ■ This makes LHC 



searches for this point difficult, as we discuss in more detail in Section IIIB 

The second point is similar to the first but has an 'intermediate-slepton' between x± an d xf ■ It 
is most like the benchmark point considered in [59]. Co-annihilation processes further deplete the 
relic density so the higgsino-fraction is smaller than before, resulting in a more massive xf an d a 
higher fine-tuning between m\ and 2m^o. We again find that the spin-dependent cross-section is 
relatively large and within reach of IceCube and XENONlT's projected sensitivity [101]. Inter- 
estingly, the spin-independent cross-section is small relative to the spin-dependent cross-section. 
The intermediate mass sleptons allow the two-body decays of X2 an d xf to open up, resulting in 



promising LHC signatures (discussed in more detail in Section IIIB). 

The third point is a A-SUSY style point with very large couplings: A = 1.5 and k = —2.18. 
We have numerically solved the 1-loop RGEs for the A, k system (including the effect of the top 
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Parameter ^ Wel1 " J Intermediate- A _ SUSY 

lempered Slepton 



A 


-0.7 


-0.7 


-1.5 


K 


-0.863 


-0.77 


-2.19 


tan /? 


4.0 


4.0 


5.45 


A x [GeV] 


-369.9 


-378.0 


-478.3 


A K [GeV] 


75.5 


74.95 


-55.9 


Mcff [GeV] 


-150.0 


-190.0 


-168.0 


Mi [GeV] 


135.0 


135.0 


128.4 


m^o [GeV] 


130.0 


133.7 


129.9 




-0.89, 0.1 


0.96, -0.06 


0.975, -0.083 


N 13 , N u 


0.39, 0.19 


-0.26, -0.09 


-0.21, 0.012 


m A [GeV] 


259.45 


267.27 


259.33 


Pl3 


0.9999 


0.9999 


0.9999 


S x = \rriA — 2m^o [GeV] 


0.55 


0.13 


0.47 


m xt ' m xt t GeV ] 


-155.8, 727 


196.8, -727 


172.8, -727 




0.999, -0.988 


0.999, -0.988 


0.999, -0.988 


m e R ^ R [GeV] 


900 


144.7 


901 


(ctu) 77 x 1CT 27 cm 3 s _1 


1.2 


1.1 


0.9 




5.0 


0.5 


0.4 




0.10 


0.12 


0.11 


cjgj x 10~ 9 pb 


1.4 


0.23 


3.1 


crg D X 10~ 4 pb 


5.4 


1.4 


0.7 


°SD x 1Cr4 P b 


4.2 


1.1 


0.5 




0.64 


0.52 


0.67 


Aa^ x 10 10 


-1.0 ±2.9 


0.8 ±2.8 


-1.4 ±2.8 



TABLE I: Properties and parameters of the benchmark points discussed in the text. All parameters have 
been defined in the text except P13, the singlet component of As- For all points, 125 GeV < mf ll < 126 GeV. 
We fix M 2 = 700 GeV and M 3 = 1.3 TcV. Other soft parameters (at SUSY scale of 1 TeV): A t = -2.2 TeV, 
A b = -1.0 TeV, A T = A^ = 0, m ic ^ = Dt = 1.0 TeV, ^ = 0.8 TeV, c (j uc b d = L6 TeV > 
;//,, = 1.2 TeV. 



quark), and find that a Landau pole occurs at around 70 TeV. Counter to naive expectations, the 
fine-tuning between tua and 2m^ is not significantly reduced. This is because A(AiVi3Afi4 — kN% 5 ) 



(c.f. Eq. (j8j)) is comparable to the value in the 'well-tempered' point. The large value of A does 
increase the spin-independent cross-section through the XSH U ■ terms (c.f. Eq. (12)) and is 
just consistent with the current XENON100 limit. The lower higgsino-fraction compared to the 
previous points results in a lower spin-dependent cross-section, although still just within reach of 
IceCube's sensitivity. 

It is well known that the Standard Model prediction for is more than 3cr below the ex- 
perimental value |102j . All of our benchmark points give a small additional contribution to the 
Standard Model value of Aa^ but not enough to significantly bring the theoretical value into better 
agreement with the experimental result. The three benchmarks points have fj, e g < while Mi and 
M2 are positive. In this case, Ref. |93] has shown that the SUSY contribution to (g — 2)« can be 
increased if the bino-higgsino-/i/j contribution dominates and there is a large hierarchy between 
the left- and right-handed sleptons. An alternative possibility for achieving a consistent value is 
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FIG. 4: The ratio of the 7Z to 77 annihilation cross-sections near the 'well-tempered' benchmark point as 
M 2 is varied. We have also varied fj, e g to keep the relic density constant and m%a = 130 GeV. The solid red 
line is the result of our analytic formula (Eq. 18) and the dotted blue line is the result of the full numerical 

1 /2 

calculation. The ~ 10% difference is due to the field renormalisation 6Z z , missing from the analytic result. 
The ratio increases as M2 is lowered, a result of the increasing wino component in the lighter chargino. 



through the introduction of flavour violation [103] but we do not consider this further. 



A. A second line from jZ 

There is strong evidence for a 7-ray line at ~ 130 GeV and intriguingly, there is also weaker 
evidence for a second line at ~ 111 GeV [U I15j . The best fit to the relative annihilation cross- 
section is (crv)nfZ I '(o"u) 77 = 0.661^43 [IB] . As we have already noted, this model does predict a 
second line with energy 



??? 



through the process shown in the right panel of Fig. ITl For m^.o = 130 GeV, the second line 
has Ery = 114 GeV. In Table [I] we quote the values ioi\av)^z/(o'v)^ n for our benchmark points, 

finding that they are all consistent with the measured value. The numbers for {av)~ f z quoted in the 

1/2 

table were calculated numerically in SloopS and include the field renormalisation 8Z z , important 
when the higgsino- fraction is large [5l] • This contribution is missing in |59} 1104] . Further technical 
details are given in Appendix [A] 

In Appendix [B] we also present an analytic calculation of (av)jz, adapting the MSSM result 
from [105J. Taking the ratio of our analytic results for (<tu) 77 and {av) 1 z, we find 



(<rv)>yz _ 1 L _ m 2 z I 5 
M 77 8sin 2 20 w { 4m 2 J Q\ ± 



~ (18) 



X? 

where J-^± and Q^± are functions that depend on the chargino masses and mixing matrices: 
m ~+ , m ~+ , U, V . The full expressions can be found in Appendix [B] From this formula, we see 



that once m^o is fixed (in our case, through the position of the 77 fine) and since mz and sin 9w 
are known, the ratio of (<tu) 7 ^/((T?;) 77 directly probes properties of the chargino sector. 
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To demonstrate this, in Fig. [4] we take our 'well-tempered' benchmark point and show how the 
ratio (crv)jz I \ av )~a changes as we vary the parameter M2. As M2 decreases, we increase /i e g in 
order that m^o and $7dm^ 2 remain fixed. We also adjust A\ and A K so that As remains singlet-like 
with mass wia ~ 2m^o. The solid (red) line is the result of our analytic formula Eq. (18) and the 



dotted (blue) line is the result of the full numerical calculation. The ~ 10% difference is due to 

1/2 

the field renormalisation 5Z z , missing from the analytic result. 

Smaller values of M2 lead to an increase in the wino-fraction of the lighter chargino and a 
smaller mass for the heavier chargino. The chargino-chargino-Z coupling is larger for a wino-like 
chargino compared with a higgsino-like chargino. Therefore, as M2 is lowered, this coupling grows 
for the lighter chargino while the chargino-chargino-7 coupling remains constant (since it is simply 
the electric charge) resulting in an increase in the ratio {crv) 1 z I '(crv) 77 . 



B. LHC signatures 

Although the singlet-like pseudoscalar Higgs with txia ~ 2m^o is essential, it does not give rise 
to any observable collider signature. Therefore, the best prospects for observation of new physics 
at the LHC in this scenario are in searches sensitive to neutralino and chargino production. All 
of our benchmark points have a relatively large higgsino- fraction in the LSP Xi- This implies that 
there will be light x\ 3 an< ^ Xii which are mostly higgsino. If there are no intermediate sleptons or 
squarks in the spectrum, these will decay via three-body decays. On the other hand, when there 
is a slepton whose mass lies between that of the LSP and the higgsinos, the higgsinos will decay 
to the slepton and then to the LSP via two-body decays. In both these cases the most promising 
signatures involve dileptons and trileptons plus missing energy carried away by the LSP. 

Both ATLAS |106|. 1107] and CMS [108J have recently published searches for direct gaugino 
production at ^/s = 7 TeV with these final states. We have implemented the ATLAS two- and 
three-lepton searches in the RIVET framework [109] (we expect the implications of the CMS results 
to be similar). We generate events using Heruig++ 2.5.2 [1 10|, 1111] and implement the experimental 
event selection using RIVET 1.8.1 to analyse the fully showered final state. The ATLAS dilepton 
search |107| does not claim sensitivity to spectra where the gauginos decay via three-body decay 
without the presence of an intermediate slepton. Our results are in agreement with this. We find 
that the ATLAS search has a low sensitivity to our 'well-tempered' benchmark point due to the 
small mass splitting (16 GeV) between the LSP and the higgsino-like charginos and neutralinos 
X2 3- This small splitting means that the transverse momentum px of the leptons is generally low, 
leading to a small acceptance in all search channels. Our benchmarks are also not ruled out by 
the trilepton searches. While we were finishing this article new gaugino searches from ATLAS and 
CMS based on data taken at 8 TeV were announced at the HCP2012 |112] 1113] conference. These 
results extend the reach with respect to m^± but do not close the gap at low m^± — m^o, which 
may require data at 14 TeV or a targeted analysis. 

To highlight the phenomenological differences between the 'well-tempered' and 'intermediate- 
slepton' cases, in Fig. [5] we show the missing energy distribution (left panel) and pr distribution 
of the the leading electron for the 'well-tempered' point (in red, solid line), and the 'intermediate- 
slepton' point (in dark blue, dashed line) at an energy of 14 TeV. We have imposed the trigger cuts 
detailed in |1U7| . along with lepton isolation and associated quality criteria. Jets are reconstructed 
using the anti-fci algorithm [114], with R = 0.4, and we only accept points that yield exactly three 
leptons. As both points have a neutralino with the same mass, the shape of the MET distribution is 
similar. However the cross-section for the well-tempered case is lower due to the smaller probability 
that three lepton pass the pt cuts. In the right-hand panel of Fig. [5] we show the pt spectrum for 
the leading lepton when it is an electron; the plot is similar for the muon case. The spectrum for the 
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FIG. 5: Left panel: The missing energy distribution for our 'well-tempered' and 'intermediate-slepton' 
benchmark points for trilepton events at 14 TeV. The distributions have a similar shape as the neutralino 
mass is the same, but the well-tempered one has a smaller cross-section. Right panel: The px distribution 
of the leading lepton (chosen to be an electron) for the same event selection. The more compressed 'well- 
tempered' spectrum is much softer than the 'intermediate-slepton' case. A targeted search with lower 
leptonic px triggers may be helpful in this case. 



'intermediate-slepton' point peaks at higher px and has a much longer tail than the 'well-tempered' 
point, as there is more phase space for the leptons in this case. Probing the 'well-tempered' region 
will prove difficult at the LHC, and may require a dedicated search with lower leptonic px triggers 
than are currently implemented. One possibility for constraining this scenario might be via a 
multilepton search, where the pair-produced gaugino system recoils off an ISR-jet, leading to three 
collimated leptons recoiling against a jet. A second possibility is to look at a dilepton plus 7 search 
since the decay X2 ~ * X?7 can be- large when X2 decays otherwise only by three-body decay |115| . 
For instance, for the 'well-tempered' benchmark point, this branching ratio is 47%. We do not 
pursue these options further at the present time. It is also worth noting that the position of the 
7-ray line in the Fermi data tells us the neutralino mass, so a very accurate reconstruction of the 
chargino and slepton masses may be possible using Mj^-based techniques |116j 1117] . for example. 

Even though it is difficult to constrain the 'well-tempered' scenario at the LHC, it is interesting 
to note that there is a complementarity between the LHC and direct detection searches for our 
benchmark points. The 'well-tempered' benchmark point is difficult to constrain at the LHC 
as the high LSP higgsino-fraction leads to a compressed spectrum. However, the large higgsino 
component of the LSP in this case will be directly probed at the next generation of direct detection 
experiments. On the other hand, when the higgsino component of the LSP is lower and the nucleon 
scattering cross-section is smaller, the splitting between the LSP and higgsinos is larger, leading to 
a harder px spectrum. Also, the intermediate slepton required to deplete the relic density leads to 
more leptons being produced in decay chains. Thus, both these scenarios are amenable to testing 
in the near future. 



IV. OTHER SINGLET EXTENSIONS OF THE MSSM 

The benchmarks points that we have discussed have all been within the NMSSM, which has 
a discrete Z3 symmetry. However, our analytic results did not rely on this symmetry so will 
apply more generally. Therefore, we briefly consider the phenomenology of a more general singlet 
extension of the MSSM. The most general extension of the MSSM by a gauge singlet superpotential 
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is 

W = Wnmssm + vHu -H d + £S+ ^ S S 2 . (19) 

One strong motivation for the NMSSM is that it solves the /x-problem of the MSSM so including 
an explicit //-term in the superpotential naively ruins this solution. However, the superpotential 
with £ = (the linear term in S can be removed by a shift in its vev) and an underlying Z^ 
symmetry leads to \i ~ \i$ ~ 0(m 3 /2) when the Z^ symmetry is broken to the usual matter parity 
after SUSY breaking. 8 This has been coined the GNMSSM (see [ml [[22] for further details). We 
briefly consider the effect of these additional terms. 

In Section [TTJ, we argued that the singlino component of the neutralino in the NMSSM is small 
for m^o ~ 130 GeV. In the GNMSSM, there is an extra contribution to ^ e ff and rrig from \i and 
\jls respectively. With these extra parameters, the singlino mass is 

m§ = 2ns + fj, s = 2 /x e ff + (ms - 2 ^m) • (20) 

Comparing with Eq. ([6]), we see that we now have the freedom to choose \i and /i<; such that 
m s ~ /"cff an d therefore, have a neutralino with m^o ~ 130 GeV where the singlino component is 
large. 

For a singlino-higgsino neutralino, we can estimate the size of the singlino component required to 
obtain the correct relic density from Fig. [3} This is because, for heavy squarks and sleptons, the bino 
component acts like a singlino component relative to the much more efficient higgsino component. 
From Fig[3j we take iV"i3 ~ 0.3, iVi4 ~ 0.1 and Ni§ ~ 0.95 in order to estimate the expected fine 
tuning required between m^o and tua- In Eq. ([8]), we considered (AA^A^— kN^ 5 ) = 0.05 and found 
5 X ~ 1.5 GeV. Using A ~ k ~ 0.7 and the numbers above, we find that (AiVi3iVi4 — kN 2 5 ) ~ 0.6, 
larger by a factor of ~ 10. As a result, we expect that S x is about an order of magnitude larger 
and the resulting tuning between m^o and rriA should be less severe in the GNMSSM. As this work 
was nearing completion, a dedicated study of the GNMSSM scenario confirms these findings |104| . 



V. CONCLUSIONS 

Motivated by recent claims of lines in the Fermi gamma-ray spectrum, we have investigated en- 
hancing neutralino annihilations into photons in singlet extensions of the MSSM. The enhancement 
occurs when the lightest neutralino is on resonance with a singlet-like pseudoscalar As (see Fig.JT]). 
Throughout, we have adopted an analytic approach as well as a numerical approach to understand 
the underlying physics. Using the analytic approach, it is clear that the lightest neutralino must 
have a large higgsino or singlino component in order to couple efficiently to this resonance through 
the XSH U • Hd or kS 3 /3 superpotential terms. 

We focused particularly on the phenomenology in the NMSSM, presenting benchmark points in 
Table|l]that are consistent with all constraints, including the continuum flux of photons and the relic 
density. In Section [TTJ we argued that the neutralino singlino-fraction is low in the NMSSM. The 
higgsino-fraction is bounded from above by the requirement of achieving the correct relic density 
but can still be as high as 25% (see Fig. [3]). A large higgsino component is correlated with a large 
nuclear scattering cross-section; both the spin-independent and spin-dependent cross-sections are 
large. In order to ensure that the spin-independent cross-section is below the current limit set by 



This J?-symmetry also cures the domain wall problem of the NMSSM [U8HT20]. 
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XENON100, it is necessary to exploit cancellations between contributions from different CP-even 
Higgs bosons. This occurs when /i e g < 0. 

In this scenario a 7Z line at energy ~ 114 GeV accompanies the 77 line. The relative strength 
of these lines depends on the chargino masses and mixings (see Fig. [4]). Currently the relative 
strength is poorly known but a more precise determination will help to constrain properties of the 
chargino sector. 

Due to the large higgsino-fraction of the lightest neutralino, in general there is a light higgsino- 
like chargino and two other higgsino-like neutralinos close in mass to the lightest neutralino. Decays 
of these particles give rise to LHC signatures with dileptons and trileptons plus missing energy. 
We investigated these signatures for two of our benchmark points; one with a slepton with a mass 
between the lightest chargino and lightest neutralino ('intermediate-slepton') and one without 
('well-tempered') (see Fig. [5]). With an intermediate slepton, the LHC signatures are promising. 
The case without an intermediate slepton is more difficult to probe at the LHC and may require a 
dedicated search with lower leptonic pt triggers than are currently implemented. 

If a 130 GeV neutralino is the source of the gamma-ray features observed by Fermi, the associ- 
ated phenomenology at colliders and direct detection experiments is potentially very rich. 
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Appendix A: Numerical implementation of (erv) 77 and (av) 1 z 

In the NMSSM, (crv) 77 and (ov) 1 z were computed in [571 [SS] by adapting the formulas of [fi9"l 
105J to the NMSSM case. A complete numerical calculation for these final states has also been 
performed in [581 1123] . The calculation of the 7Z final state requires the calculation of the field 
renormalisation 5Z^ which is generated from the (tree-level) XiXi^ vertex through a Z — 7 
one- loop transition. This contribution was missing in [105] (and by extension [59]) and is gauge- 
dependent, as was first pointed out in [H] for the MSSM calculation of (av)~ / z- Besides being 
needed to obtain a gauge invariant result, this contribution can be numerically significant when the 
neutralino has a significant higgsino-fraction, as is the case we study here. This can be understood 
by the fact that the coupling of the neutralino to the Z boson is proportional to the neutralino's 
higgsino component. 

Our implementation of the cross-section (<ri;} 77 follows [58], which relies on the SloopS code |54^ 
[7TH73] . The SloopS code benefits from a non-linear gauge- fixing which enables the user to check 
the gauge invariance of the result. Further details about the numerical implementation for the 
NMSSM calculation can be found in |58| . Note that by default no width is implemented in SloopS 
since it can spoil gauge invariance. As the width of the singlet As is very narrow the inclusion of 
the width is completely negligible for the points we consider. Although the inclusion of the width 
in principle relevant when very close to the resonance 2m^o = rriA, these points give a rate that is 
much too large and would already be excluded, as discussed, for instance, in [58]. In addition to 
the diagram of Fig. [l] there are contributions where the charginos running in the loops are replaced 
by leptons and quarks. Analytically their contribution is similar to the formulas in Eq. ^ where 
the term AL42V12 is replaced by the Yukawa coupling of the corresponding lepton/quark. However, 
unlike the chargino, these diagrams are suppressed since mi >q 56 m^o. These diagrams are therefore 
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sub-dominant. There is also a slight destructive interference between the resonant type diagrams 
and box-type ones with charginos inside. However, these are also numerically sub-dominant. 



Appendix B: Analytic expressions for (av) x °x1^A s ^-yy an< ^ ( av ) 'xlxl^A s ^-yZ 

In this appendix we present the analytic expressions for (<ru) 77 and (av)~ f z- : considering only 
the contribution from the s-channel psuedoscalar resonance As, which we assume is a pure singlet 
state. To derive these results, we adapted the calculations from [69] and |105j for (av ) 77 and {<rv)yz 
respectively. Unlike our numerical result, our analytic result for (av) y z does not include the field 
renormalisation 5Z 1 ^ . 

-yZ 

Including both charginos in the triangle and assuming that A and k are real, we find that 

a 2 \ 2 (\N 13 N 14 - kN 2 5 ) 2 n2 

Wnt ~ 47r 3 (4m 2 o _ m 2 )2 + T 2 2 y& > W 

where we have defined 



G%± = ^ rn~+Uj2Vi2 arctan 2 



2 

m± 



and U and V are the usual unitary matrices that diagonalise the chargino mass matrix. 
Similarly, we find 



2 1 A m 2 z \ \ 2 {\N 13 N m -kN? 5 ) 2 



Here, 



where 



(B2) 



{(Jv) ~ iZ 32vr3 sin 2 26 w I 1 4m 2 J (4m 2 - m^) 2 + T\m\ ' (B3) 

\ zK.\ / Xi 

where we have defined 

~ 1. m^SijI x -2-, -2-, 1, — 2- + m^+Dijl-t ,1 . (B4) 



mi+ mz+ 4m_ . I x i \ m_ + m_ + 4m 

V V V / \ V V V 



^ = ^2^2^ + UflVvdj? (B5) 

Ai = UaVfrdjf + C/ i2 ^ 2 OJ (B6) 



= -VjlVfi ~ \v j2 V& + Sij sin 2 ^ (B7) 
6 if = -UjiUil-\ U-2 Ui2 + <% sin 2 9 W . (B8) 

Finally, we also have 

f 1 dx 

lf(a,b,c,d)= / — [slog(— 4a, 6, c; x) — slog(— 4d, b, c; x)] , (B9) 
Jo x 

where 

slog(a, 6, c; x) = log [| — ax 2 + (a + 6 — c)x + c|] . (BIO) 
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Once m^o is fixed (and since mz and sin 6w are known) the ratio 



2 



(cri;) 77 8sin 2 20 w \ 4mi / Qi 



(Bll) 



depends only on the chargino masses and mixing matrices through Q^± and J~y±- Therefore, an 
accurate measurement of this ratio will constrain properties of the chargino sector. 



Appendix C: Direct Detection Results 

Here we briefly review the computation to calculate the neutralino scattering rate at direct 
detection experiments. The low energy neutralino-quark effective Lagrangian is 

£eff = aqXiXiQQ + dgxh^^XiQl^q , (CI) 

where we only consider contributions that are not velocity or momentum transfer suppressed. 
The first-term gives rises to a spin-independent interaction while the second gives rise to a spin- 
dependent interaction. We consider each in turn. 

For spin-independent interactions, experiments typically quote the cross-section to scatter off a 
nucleon, given by 

4 ^o {Zfp + { A - Z)f n f 
= — & • (° 2 ) 

where /Ujv^o is the neutralino-nucleon reduced mass and A and Z are the nucleon number and 
charge of the target nucleus. The couplings to protons and neutrons, f p and f n , are related to the 
quantities entering C e s by 



fN r, 

— -2^ —H + w'Ta — ■ W 

q=u,d,s q Q=c,b,t v 



Here, the fx are determined from lattice QCD and chiral perturbation theory and 
fr G = 1 — J2q=uds fc q |124j . Unless stated otherwise, we take cto = 35 MeV and ovjv = 45 MeV, 
leading to 

f^ u = 0.020, = 0.026, = 0.13 (C4) 

/?„= 0.014, ff d = 0.036, #, = 0.13 (C5) 

These values are comparable to the values used in DarkSUSY |125j . The uncertainty associated 
with the strange quark matrix element fx s is large; more recent values are typically much lower. 
For instance, Ref. |126| find fj, = 0.053 which leads to a cross-section that is smaller by a factor 
of a few (comparable values of /t s were also found in [127H130] ). 

In principle, a q receives contributions from the t-channel exchange of CP-even Higgs bosons and 
s-channel squark exchange. However, for the discussion below, we ignore the contribution from 
the squarks because of the strong LHC limits on first two generations of squarks |131| I132j (in the 
absence of unusual features in the spectrum such as compression or i?-parity violation which are 
not present in our model). In this case, for u-type quarks, 

^= 9 o y 9 Jv°&i. (C6) 

m u Amy/ sin p ^— ' mi 

j=l n i 
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The corresponding expression for ad/ ma can be found by replacing sin/3 — > cos/3 and Sj 2 — > Sj±. 
The matrix Sij diagonalises the CP Higgs mass matrix such that the mass eigenstates hi, ordered 
such that m/ ll < < rrih 3 , are related to the weak eigenstates hj cak = (HdR, H u r, Sr), where 
\f2H a = H a R + iH a j, through hi = SijhJ eak . The coupling gh tX x f° r * ne lightest neutralino is 

g hiXX =g(N 12 - tan^JVn)(5iiJVi3 - S i2 N u ) (C7) 
- V2X(S a N 14 N 15 + S i2 N 13 N 15 + S i3 N 13 N u ) (C8) 
+ V2 K S i3 Nf 5 . (C9) 

The first line reduces to the MSSM contribution while the second and third lines come from the 
XSH u Hd and |5 2 superpotential terms respectively. 

We next consider the spin-dependent interaction. The additional terms do not give rise to any 
new terms not suppressed by velocity or the momentum-transfer so the result is the same as for 
the MSSM. The dominant contribution comes from t-channel Z-exchange so couples directly the 
higgsino component. The cross-section is 

12/i 2 -o 

SD _ ^P:"Xl 




<',: i w i . (cio) 

where 

d « = W*? - l*4| 2 ) (CID 

and T 3q is the third component of hypercharge. Aq' n measure the fraction of the nucleon spin 
carried by a given quark and we use 

A p u = A3 = 0.842 = A™ = -0.427 A§ = A™ = -0.085 (C12) 

With these numbers, we find that the cross-section to scatter off a proton is larger than the 
cross-section to scatter off a neutron 

<xf « 4.0 x 10- 4 pb ^I3l 2 -I^i4l 2 ^ 2 (C13) 
o*> - 3.1 x 10- pb ( |JVl3|2 7 1 |jVl4 ' 2 ) 2 • (C14) 
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